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Abstract 

The hot electron-positron unmagnetized plasma with small fraction of cold 
electron- ion plasma is investigated. The modulational interactions of the electro- 
magnetic waves and the electron-sound waves are studied. The possibility of soliton 
formation is investigated. 



I. Introduction 

During the past few years considereble amount of papers have been devoted to the 
analysis of electromagnetic (EM) wave propagation in hot, pure electron-positron (e-p) 
plasmas since e-p pairs are thought to be a major constituent of the plasma emanating 
both from the pulsars and from the inner region of the accretion disks surrounding the 
central black holes in active galactic nuclei (AGN) (Michel, 1982; Begelman et al, 1984). 
Such a plasma is formed also in the early universe (Rees, 1983; Tajima and Taniuti, 
1980). Although the relativistically hot e-p pairs form most of astrophysical and cosmic 
plasmas, a minority of cold electrons and heavy ions is likely to be present (Berezhiani and 
Mahajan, 1995). For instance, outflows of e-p plamsa from pulsars entering an interstellar 
cold, low density electron-ion (e-i) plamsa forms two temperature electron-positron-ion 
(e-p-i) plasma. The three-component plasmas - hot e-p plasma with small fraction of 
heavy ions - have been studied in the context of pulsar magnetospheres by Lakhina and 
Buti (1981) and by Lominadze et al (1986). The creation of stable localized structures of 
relativistically strong EM radiation in hot e-p-i plasma have been shown by Berezhiani 
and Mahajan (1994,1995). 

In two temperature e-p-i plasma the interesting phenomena differing from that of one 
temperature can exist. It is now believed that strong monochromatic waves emitted by 
pulsars are subject to parametric instabilities even in quite underdence plasmas. In this 
context in present paper we consider the propagation of strong EM radiation in a hot e-p 
unmagnetized plasma with small fraction of cold e-i plasma. We show, that the presence 
of a minority of cold electrons and ions can lead to the scattering of the pump EM wave 
into the electron-sound and EM wave; to the instability of hot e-p plasma against the low 
frequency (LF) perturbations. Hence, in contrast to the case of the pure e-p plasma, in two 
temperuture e-p-i plasma the three wave decay instability may occur. The possibility of 
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the soliton formation due to the modulational instability of EM wave is also investigated. 

II. The model and formulations 

Let us assume that the velocity distribution of particles is locally a relativistic Maxwellian. 
Then the dynamics of the fluid of species a (a — e,p,i) is contained in the equation 
(Javakhishvili and Tsintsadze, 1973): 

■^-(KU ak W a ) - -^-P a = -F lk J ak , (1) 

where U l a = [j a , r y a u a /c] is the hydrodynamic four velocity, u a is the hydrodynamic 
three- velocity of the hot e-p fluid, 7« = (1 — u 2 a /c 2 )~ l l 2 is the relativistic factor, J ak is 
the four current, F lk is the electromagnetic field tensor and W a is the enthalphy per unit 
volume: W a = (n a /j a )m a c 2 G a [m a c 2 /T a }. Here m a and T a are the particle rest mass 
and temperature of species a, respectively, n a is the density in the laboratory frame of 
the e-p-i fluid and G a (z a ) = K 3 (z a )/ K 2 (z a ), (z a = m a c 2 /T a ), where K v are the modified 
Bessel functions. For the nonrelativistic temperatures (T a « m a c 2 ) G a = l+5Ta/2m a c 2 
and for the ultrarelativistic temperatures (T a >> m a c 2 ) G a = 4T a /m a c 2 >> 1. The 
relativistic pressure in the rest frame is P a = {n a /^ a )T a . 

We assume that hot electron and positron temperatures are equal and constant while 
the process of EM wave interaction with given fluid (G a h = const). Note that here and 
below the subscript "c" is used for cold electrons and "h" - for hot particles respectively. 

^From the set of equations (1) the equation of motion can be written as follows: 

^(P a G a ) + —VP a = e a E + —(u a x B), (2) 

at n a c 

where P a = r y a m a u a is the hydrodynamic momentum, E and B are the electric and 
magnetic fields and d a /dt = d/dt + u a ■ V is the comoving derivative. For cold electrons 
in the eq. (2) G c — 1 and T ec = T C = const should be assumed. 
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And for all kinds of species we have the continuity equation: 



y + VM = 0. (3) 

To study the nonlinear propagation of intense EM wave in a relativistically hot e-p 
plasma with small fraction of cold e-i plasma we must couple the eq.-s of motion with the 
Maxwell equations. In the terms of the potentials defined by: 

E = --— -V0; B = VxA, (4) 
they take the form (Coulomb gauge V • A = 0): 

- c 2 AA + c^(V0) - 4ttcJ = 0, (5) 

and 

A0 = -4np, (6) 
where for the charge and current densities we have respectively: 

P = J2 e * n a\ J = e a n a u a . (7) 

a a 

The equilibrium state for hot e-p plasma is characterized with charge neutrality (with 
unperturbed number densities of the hot electrons and positrons equal to n oh ). For small 
fraction of cold e-i plasma the equilibrium state is characterized also by charge neutrality 
(with background ion density n oc ) and 

n oc « n oh . (8) 

Also we assume that: ions are immobile; in equilibrium state hot electrons and 
positrons have the same temperatures equal to T Q h and 
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Toh » T c ; 



0. 



(9) 



Let us analyse the one-dimensional propagation (J^ 7^ 0, ^ = 0, = 0) of circularly 
polarized EM wave with a mean frequency w Q and a mean wave number k a along the z 
axis. Thus 



where A(z,t) is a slowly verying function of z and t and x and y are the standard unit 
vectors. The gauge condition gives us A z = 0. Then the transverse component of eq.-s of 
motion (2) are integerated yielding: 



where the constant of integration is set equal to zero since particle hydrodynamic moments 
are assumed to be zero at the infinity where the field vanishes. 

Now it is necessary to write the equations for longitudinal motion. This motion is 
driven by the ponderomotive pressure (~ P£ ± ) of high frequency (HF) EM fields and 
latter doesn't depend on the particle charge sign. In purely e-p plasma since the effective 
mass of the electrons and positrons are equal (G e — G p — G) the radiation pressure gives 
equal longitudinal moments to both the electrons and positrons and effects concentration 
without producing the charge separation (n e = n p and = 0) (Berezhiani and Mahajan, 
1994; Kartal et al, 1995). But as it was shown by Berezhiani and Mahajan, (1995) the 
introduction of small fraction of heavy ions leads to "symmetry breaking" between hot 
electrons and positrons and it becomes possible to have finite 0. As we will see below, 
due to the presence of small fraction of cold e-i plasma in hot e-p plasma the electrostatic 
potential will be surely created. 




(10) 




(11) 
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Let us redefine the electron rest mass in equations for hot e-p plasma as: 



m^mG h (T h ) = M.. (12) 
Then the eq.(5) for transverse motion can be written as: 

d 2 2 d 2 2 , A ± 2 n c A_l 

— A ± - c — A ± + u h (n eh + n ph ) h uj c = 0, (13) 

dt dz 2 n oh >y h n oc 7 C 

where uj\ = 4iie 2 n h/M, oj 2 = 4ne 2 n oc /m e and 7^ and 7 C are the relativistic factors for 
hot and cold electrons respectively, 7 Q = (1 + P^/m^c 2 ) 1 / 2 ; n c , n e h and n p h are the 
cold and hot electron and positron densities respectively. 
The equations for longitudinal motion have the form: 

dt dz \m a >y a J 

(- + Pza d ) P za + ~- (-) = -e a ^ - Cl dlAl \ (15) 

where for hot particles m a h = M is assumed and in P za h the mass redefinition is performed 
(see the relation (12)). 

In what follows we consider only the weak relativistic case assuming P 2 /m 2 c 2 << 1. 
In the presence of two different temperature electron plasma for LF motion it is possible 
to satisfy the condition (Khirseli and Tsintsadze, 1980): KV Tc « Q « KV Th , where 
Vtc and Vtk are the cold and hot electron thermal velocities respectively and Q" 1 and 
K^ 1 are the characteristic time and spacial spreads of the pulse (Q « uj q ; K « k Q ). 

First let's find the equation for LF motion. Under the above mentioned conditions for 
hot particles we have: 

c d<j) _ T h d (n ah \ | e 2 8\A\ 2 
a dz n ah dz\j h ) 2M-f h c 2 dz 
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From eq.-s (16), introducing 5n c = n c — n oc and 5n a h = n a h — n h (Sn c « 
; 5n ah « n oh ), we obtain: 

e 2 \A\ 2 e 2 \A\ 2 
Sn h = 5n ph + 5n eh = -n oh j^-^ + n oh ^^j. (17) 

Using equations (6), (7) and (17) we find the relation between 5n e h and 5n c written 

as: 

]_ e 2 \A\ 2 e 2 \A\ 2 

8n e h = — 5n c — n h——— — - + n ft— -— — . (18) 
2 2MT h c 2 2M 2 c 4 v ; 

Using the eq.(18) in the eq.(16) written for electrons finally we obtain: 

d± = _Jj^d_ 5nc (1Q) 

dz 2en oh dz 

Thus, as we already mentioned above, the presence of small fraction of cold e-i plasma 
in hot e-p plasma gives rise to electrostatic potential. 

Substituting the eq.(19) in the eq.(15) written for cold electrons, after simple algebra, 
assuming that: 



^ « - — , (20) 
T h 2 n oh 

one can get the equation for 5n c : 

d 2 d 2 e 2 d 2 

—Sn c - c 2 —5n c = n oc ^ 2 ^\A\ 2 , (21) 

where 

' 1 n oc T h \ 1 I 2 



C 2 s 



-) 



,2n oh m. 

is the so-called "electron-sound" velocity. Thus, due to the fact that the most part 
of electrons are relativistically hot, and consequently heavy (Gh ^ 1) than the small 
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part of cold electrons, it is possible to induce the "electron-sound" wave; the exciting 
ponderomotive force is defined by the HF pressure on cold electrons. 

The eq.(21) is coupled with the following equation for A (for HF wave frequency 
satisfying the dispersion relation: uo 2 = k^c 2 + 2ujI + a>g): 



The system of equations (21), (22) together with the relation (17) describes the non- 
linear wave dynamics in a relativistically hot e-p plasma with small fraction of cold e-i 
plasma. As we see we have the scattering of EM pump wave into the electron-sound 
and EM wave. Note that in the purely e-p plasma the three wave scattering processes 
do not exist. The presence of small fraction of cold e-i plasma here is a reason to have 
the LF longitudinal waves together with the HF EM waves. Using this result we can 
conclude that the radiation emanating both from the pulsars and AGN entering the cold 
low density e-i plasma undergoes the modification due to the scattering processes. 

It is possible to find the stationary solution of the system of eq.-s (17), (21), (22). We 
look for the solutions as: 





(22) 



where v g is the group velocity of HF wave. 



A = A(£,t); 5n c = 5n c (£,r); £ = z - v g t; t = r; — << v g — . 



(23) 



In the subsonic regime: v g << c s , i.e. 




(24) 



from the eq.(21) we obtain: 
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and substituting the eq.(25) in the eq.(22), if we have the relativistically hot e-p plasma 
(G h » 1), we get the Nonlinear Shrodinger Equation (NLSE): 




■A = 0. 



(26) 



As it is wellknown the eq.(26) has the stationary solution representing the subsonic 
soliton of rarification (the total density variation 5n = 5n h + 5n c < 0). 

In the case of the nonrelativistically hot e-p plasma (Gh — 1) and k Q — > we aslo 
get the soliton solution of obtained NLSE. Under definite conditions it is possible to have 
the supersonic solitons too. It should be noted that for the EM waves with v g = in 
pure, hot e-p plasma the possibility to have the stable soliton-like structures was found 
by Kartal et al, 1995. 

Let's investigate the stability of two-temperature e-p-i plasma. For this reason we look 
for A and 5n c as: 



where a(z,t) and 6(z,t) are the slowly varying in time and space real functions and 
5a « a a , 56 « 9 . 

Linearizing the system of equations (21), (22) using the relation (17) we easily obtain 
the dispersion relation: 



A(z, i) = a(z, t)e %e ^ z '^; 5n c = 5n c exp[ikz — iut] + c.c. 



a = a Q + 5aexp[ikz — iut] + c.c; 9 = O + 59exp[ikz — iut] + c.c, 



(27) 
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+ (28) 

which in the limit M » m (i.e. Gh » 1, the relativistically hot e-p plasma) for the 
coinciding roots: 

v' k 2 

u~c s k + iT~ v g k + + iT (29) 
gives the relation for increment: 



r^c,«. (30) 

Thus, the addition of even very small amount of cold e-i plasma (n oc ^ 0, i.e. c s ^ 0) 
leads to the instability of hot e-p plasma against the LF perturbations. Such three wave 
decay instability doesn't exist in pure e-p plasma. The present result shoud be useful to 
understand the character of the pulsar and AGN radiation. 

In conclusion, we have shown that in the hot e-p plasma with small fraction of cold 
e-i plasma, it is possible to have the scattering of EM wave with relativistically strong 
amplitude into the longitudinal so-called "electron- sound" and EM wave. Under the 
definite conditions the possibility of soliton solution creation for EM wave is found. 
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